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We have constructed a spherically symmetric structure model in a cosmological background filled 
with perfect fluid with non-vanishing pressure as an exact solution of Einstein equations using the 
Lemaitre solution. To study its local and quasi-local characteristics including the novel features of 
its central black hole, we have suggested an algorithm to integrate the equations numerically. The 
result shows intriguing effects of the pressure inside the structure. The evolution of the central black 
hole within the FRW universe, its decoupling from the expanding parts of the model, the structure 
of its space-like apparent horizon, the limiting case of the dynamical horizon tending to a slowly 
evolving horizon, and the decreasing mass in-fall to the black hole is also studied. The cosmological 
redshift of a light emitted from the cosmological structure to an observer in the FRW background 
is also calculated. This cosmological redshift includes local and cosmic part which are explicitly 
separated. We have also formulated a modified NFW density profile for a structure to match the 
exact solution conditions. 

PACS numbers: 95.30.Sf,98.80.-k, 98.62.Js, 98.65.-r 


I. INTRODUCTION 

The term cosmological black hole (CBH) is used to describe a collapsing structure within an otherwise expanding 
universe after the radiation era. This is different from the term astrophysical black holes (ABHs) coined to use for the 
applications of static or stationary asymptotically flat black holes in astrophysics (see for example recent paper Q)- 
CBHs are fundamentally different from the ABHs in being non-stationary and having a dynamical horizon, although 
being as much interesting as ABHs for the astrophysical applications and well suited to be used as a test arena for 
fundamental physics in strong gravity and quasi-local phenomena in weak gravity regimes. Despite extensive studies 
in relativistic structure formation and static and stationary black holes we still have little information based on exact 
general relativistic studies regarding the main features of such overdense cosmological regions including its central 
black hole. 

Since the early beginning of the discovery of the expansion of the universepeople have been looking for models 
describing an overdense region in a cosmological background ([2[, see also Q). In the matter dominated era the 
initial expansion of an overdense region within an expanding universe will finally decouple from the background and 
collapses to a dynamical black hole, in contrast to the overdensity regions in the radiation dominated era in the early 
universe where the density perturbation has to be of the order of the horizon to collapse to a primordial black hole 
(PBH). In the matter dominated era, we expect a CBH within the resulting structure and a very weak gravitational 
field outside it, being different from the familiar Schwarzschild one in being neither static nor asymptotically fiat 
[1|. Therefore, such cosmological structures, if based on exact solutions of general relativity and not produced 
by a cut-and-paste technology, are very interesting laboratories to study not only general relativistic non-linear 
structures, their quasi-local features such as mass and horizons Q, black holes thermodynamics and information loss 
puzzle [il, but the validity of the weak field approximation as well Q. After all, the universe is evolving and 
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asymptotically not Minkowskian. Therefore, one needs to have a dynamical model for a black hole, to be compared 
with the familiar results in the literature on black holes within a static and asymptotically flat space-time [H, where 
global concepts such as the event horizon are not defined. Nevertheless, depending on the specific cosmological model 
we may be able to define a hypersurface separating the trapped null outgoing geodesics from those approaching 
(future-)infinity. This hypersurface effectively plays the role of the event horizon y. Numerically, this hypersurface 
is determined by tracing the outgoing null geodesics at late times; We will refer to it as the quasi event horizon. 

The need for a local definition of black holes and their horizons has led us to concepts such as Hayward’s trapping 
horizon [9[, isolated horizon [T^, Ashtekar and Krishnan’s dynamical horizon (DH) [ll|, and Booth and Fairhurst’s 
slowly evolving horizon El- The CBH we are going to study is an excellent example of testing these different 
concepts and their relationship in addition to understand its difference to an asymptotically Minkowskian and static 
black hole. Now, a widely used metric to describe the gravitational collapse of a spherically symmetric dust cloud 
is the so-called Tolman-Bondi-Lemaitre(LTB) metric [l^- Exact general relativistic models for the dynamic of an 
asymptotically FRW structure leading to a central dynamical black hole based on a LTB metric with no pressure 
has been reported in im. These models may be extended to a perfect fluid with a non zero pressure, the so-called 
Lemaitre models [I 3 . Our interest is now using these inhomogeneous cosmological solutions of Einstein equations 
as a model for a cosmic structure with non-vanishing internal pressure leading to a dynamical CBH within a ERW 
matter dominated expanding universe. We, therefore, try an ideal fluid with an interpolating pressure function 
being non-zero inside the structure and vanishing at infinity where we expect a matter dominated ERW universe. 
To achieve this goal, we have to avoid any cut-and-paste technology of finding the solution. Any internal solution 
pasted to a FRW universe does not reflect the dynamics of the inhomogeneous universe due to the homogeneity of 
FRW universe outside the structure. In contrast, our model is just asymptotically FRW reflecting the full relativistic 
local and quasi-local effects due to the cosmic fluid. The ideal fluid is modeled such that the non-vanishing pressure 
inside the structure goes smoothly to a matter dominated universe far from the structure. A similar problem studied 
extensively in literature in the last 40 years is the issue of primordial black holes (PBHs). These are structures 
within the radiation dominated phase of the universe, usually in the late phase of inflation, due to the superhorizon 
perturbations [T8|. PBHs are usually based on the same Lemaitre cosmological solutions of Einstein equations which 
are spherically symmetric and inhomogeneous. The perturbation is formalized by assuming that the horizon size Rh 
is smaller than the structure size or by assuming e = ^ < 1, where the relevant quantities are expanded in 
the powers of e. Using the same Lemaitre solution although not in the same coordinates, we are interested in cases 
where e ^ 1, i.e. structures are much smaller than the horizon size as expected for cosmological structures within 
the matter dominated FRW universe. Therefore, we expect to recover novel features and we have also to propose 
a new algorithm how to solve numerically the field equations. In addition we are asking different questions such as 
the kind of black holes we may encounter and the characteristics of their apparent and event horizons, the effect 
of the pressure inside the structure on the collapse behavior and the matter flux. There are many other questions 
encountering in black hole literature to be faced in future studies, such as the definition of a non-rotating spherically 
symmetric dynamical black hole within an expanding universe, its differences to the simple Schwarzschild or Kerr 
black hole, its event and apparent horizon and their features, the internal structure of CBH, the rate of collapse, 
the effect of pressure inside the structure and its probable non-Newtonian and non-linear novel effects due to the 
non-vanishing matter outside, and the information puzzle. That is why the study of CBHs goes beyond the study of 
PBHs and is a new arena for a novel black hole terminology. 

Section H is an introduction to the spherically symmetric inhomogeneous prefect fluid cosmological models. In 
section HI the result of the numerical integration is reported expressing the main characteristics of our model 
assuming different pressure profiles. In the section IV, the redshift of the light coming from nearby the structure is 
investigated. In section V we develop a modified NEW density profile for a cluster of galaxies taking into account 
the requirements of our solution and study its evolution. We will then discuss the result in section VI. Throughout 
the paper we assume SttG = c = 1. 


II. GENERAL SPHERICALLY SYMMETRIC SOLUTION 


Consider a general inhomogeneous spherically symmetric spacetime Ea filled with a perfect fluid and a metric 
expressed in the comoving coordinates, = (t, r, 0): 

ds^ = — dt^ + dr^ + dQ^ , 


( 1 ) 
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where a = cr{t,r)^ X = A(t,r) are functions to be determined, R = R{t,r) is the physical radius, and dft‘^ = 
d0‘^ + sin^ 0 dcj)^ is the metric of the unit 2-sphere. The energy momentum tensor of the perfect fluid is given by 

= {p^p)u^^u^ ( 2 ) 

where p = p{t^r) is the mass-energy density, p = p{t,r) is the pressure, and = (e“^, 0,0,0) is the perfect fluid 
four-velocity. 


A. Field Equations 

In addition to the Einstein field equations, A^ we will use the conservation equations in the form 


26^^ 

(p+p) 


v„r*'' = A + 


2p AR 

- - - 1 - =0 

ip + p) R 


( 3 ) 


{p + p) 


= a' 


P 


P + P 


= 0 , 


( 4 ) 


where the dot means the derivative with respect to t, and the prime means the derivative with respect to r. The 
Einstein equations lead finally to following equations 


and 


_9 

dr 


dt 


R + i?^2g-2<T _ 

3 


R + RlTe-^^ - RR'^e-^ - -AR^ 

o 


= hvpR^R' 


= —npR^R. 


The term in the brackets is related to the Misner-Sharp mass, M, defined by 


= i?2e-2<T _ i?'2g-A ^ 

R 3 


Eqs m and ([6]) may now be written as 


K,p = 


2M' 


2M 


R^R' R^R 

We may write Eq © in the form of an evolution equation of the model: 


«=Wf 


where 


f{t,r)=R'^e 


f2-X 


( 5 ) 


(6) 


( 7 ) 


(8) 


( 9 ) 


( 10 ) 


is the curvature term, or twice the total energy of test particle at r (analogous to /(r) in the LTB model). In this 
paper we have set A=0. Note that R{t, r) can not be directly obtained from this equation because of the unknown 
functions A, a, and M. 

The metric functions gu and grr niay be obtained by integrating ® and ( 0 : 


cr = c{t) - f 

Jr, 


p' dr 
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( 11 ) 


( 12 ) 
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where c(t) and Aq (r) are arbitrary functions of integration (see 0 for more details). In the case of c(t) it is easily seen 
that requiring our coordinates to lead to the LTB synchronous ones for p = 0 leads to c{t) = 0. We notice also that 
according to m and the LTB coordinate conditions, the choice of Ao(r) is equivalent to the choice of /(to, r) = fo{r). 
One may prefer to choose fo{r) and then calculate Ao(r) from = Rq /{I + /o)- 

We have therefore 5 unknowns p, p^cr, X, and four dynamical equations p, M, A, and i?, in addition to an equation 
of state p = p(p), and the definition of the mass M ([7]). This defines a numerical algorithm to find solutions for the 
dynamics of our spherical structure after assuming the initial conditions. The algorithm for solving the coupled PDEs 
can be found in 


III. EQUATION OF STATE AND THE RESULTS 

We are now ready to specify the equation of state and integrate the model to see its characteristics. To have a 
comparative discussion of the results we consider two types of equation of state: a perfect fluid with a constant state 
function, p = wp^ and a more general case with the equation of state p = ws{r)p matching our needs for a structure 
with pressure inside the structure and a pressure-less matter dominated universe far from the structure. We may 
then choose the function s{r) in a way that the pressure becomes zero at infinity, i.e. at r >> ro- A suitable choice 
is s{r) = e~^ where ro is the distance of the void (boundary of expanding and collapsing phase) from the center of 
the structure. This is a more realistic model to describe a black hole collapse within the FRW universe and to see 
the effect of the inside pressure while the universe outside is matter dominated with no pressure. 

The model we envisage starts from a small inhomogeneity within a FRW universe. The density profile should 
be such that the metric outside the structure tends to FRW independent of the time while the central overdensity 
region undergoes a collapse after some initial expansion. At the initial conditions, where the density contrast of the 
overdensity region is still too small, we may assume that the metric is almost FRW or LTB; the density contrast and 
the pressure does not play a significant role. The dynamics of Lemaitre universe will give us anyhow the expected 
structure at late times. To choose the initial conditions at the time to, we will therefore use a LTB solution with 
a negative curvature function. We have in fact tried both examples of LTB or FRW initial data and received no 
significant difference between the final Lemaitre solutions. 

Now, let us choose the the two initial functions /(r) = /(to,r) and M(r) = M(to,r) in the following way to achieve 
an asymptotically FRW final solution: 


/(^) = 


M{r) = 

Far from the central overdensity region we have 

linir^^fir) = 0 , 


(13) 

(14) 

(15) 


li/TTlj^ —( t *) ^ 


(16) 


showing the asymptotically FRW behavior of the initial conditions. The corresponding LTB solution of Einstein 
equations now gives us R{r) = i?o(?’) at the initial time to- By choosing suitable a and b the ’’reality condition” 0 


2M 


KR^ 


> 0 , 


(17) 


will be satisfied. This condition can be obtained from Eq ([9]). In addition to ’’reality condition”, the ’’weak energy 
condition” j2lj must be satisfied at every r and t; 

p(r, t) > 0, p(r, t) + p(r, t) > 0, 

It can be seen from Fig. m that our initial conditions give rise to p(r, t) > 0 at every t and r. 


( 18 ) 
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(a)Density evolution of our cosmological black hole for the perfect (b)Density evolution of our cosmological black hole for the perfect 
fluid with the equation of state p = wp. fluid with the equation of state p = wps(r). Note the less signiflcant 

central density and the more flat density profile near the center of 
the structure. 


FIG. 1: Density evolution of our cosmological black hole. 


Assuming an equation of state to satisfy weak energy condition, p{r, t) > 0, is now enough to numerically calculate 
the necessary dynamical functions of the model. Specifically, by looking at R{t^ r) and p(t, r) we may extract 
informations of how the central region starts collapsing after the initial expansion and how a black hole with distinct 
apparent and event horizons develops while the outer region expands as a familiar FRW universe. We may also find 
out the difference to the case of the pressure-less model. It will also show if and how the very weak gravity outside 
the collapsed structure affects the dynamic of the central structure in comparison to the familiar Schwarzschild 
model. The results of the numerical calculation for both equation of states are given in the following sections. 


A. The density behavior 

The density profiles for both equation of states as a function or t and r are given in Fig.([T]). A comparison of 
these figures shows the effect of the pressure on the development of the central black hole. Obviously in case of 
non-vanishing pressure outside the structure the collapse is more highlighted with a more steep density profile. The 
over-density region in the collapsing phase is always separated from the expanding under-density region through a 
void not expressible in these figures. We will consider the deepest place of the void as the boundary of the structure. 
This boundary is always near by the boundary of the contracting and the expanding region of the model structure. 

In the numerical calculation at late times the density near the center of the black hole approaches infinity leading 
to numerical errors and crashing of the algorithm. We refer to this effect as the black hole singularity, as discussed in 
HHi. The cosmological singularities (big bang and big crunch) are not included in our numerical simulation due to 
the initial conditions being fixed at a time much later than the big bang in the matter dominated era. Assuming the 
background to be the standard flat FLRW metric, there is no big crunch singularity in our model. By selecting other 
backgrounds such as the closed FLRW model, the big crunch singularity will play a crucial role in the development 
of singularities Q. 


B. The pressure effect 


Figs.(2(a)) and (3(a)) show the behavior of the collapsing and the expanding regions for the equation of state 
p = wp by depicting the corresponding Lemaitre Hubble parameter R/R versus the physical radius R. Figs. (2(b)) 
and (3(b)) show similar data for the equation of state p = ws{r)p. Note that the function s(r), as defined to get matter 
dominated FRW universe at far distances, has no significant effect, and the qualitative behavior of the dynamics of 
the physical radius is independent of it. Therefore, as far as we are interested in the qualitative features of the model, 
we will just use the simple equation of state with s{r) = 1. 

The place of separation between the expanding and collapsing region defined by R > 0 and R < 0 is almost coincident 
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(a)The overall scheme: R > 0 and R < 0 show the expanding and (b)A magnified scheme related to the regions near to the center 

collapsing regions, respectively. 


FIG. 2: The behavior of the Hubble parameter R/R in the case of p = wp. Evidently the pressure slows down the collapse 
velocity near the center of the structure 



(a)The overall scheme: R > 0 and R < 0 show the expanding and (b)A magnified scheme related to the regions near to the center, 
collapsing regions, respectively. 


FIG. 3: The behavior of the Hubble parameter R/R in the case of p = wps{r). The features are qualitatively as in Fig ([2]). 


with the place of the void where we have defined as the boundary of the structure. Now, from the figures we realize that 
the effect of the pressure in different regions of the model and its comparison to the homogeneous FRW model is an 
intriguing one. As we know already from the Friedman equations in FRW models, the pressure adds up to the density 
and has an attractive effect slowing down the expansion leading to a more negative acceleration = — |(p + 3p)). 
This is evident from the figures at distances far from the center where our model tends to an FRW one. Whereas 
within the structure where we have a contracting overdensity region the behavior is counter-intuitive. Except for the 
case of vanishing pressure, in all the other cases the pressure begins somewhere to act classically like a repulsive force 
opposing the collapse of the structure. To see this more clearly, we have also depicted the acceleration in Fig.®. 
As we approach distances near to the center, the negative acceleration in the FRW limit and even inside the void, 
gradually increases to positive values, meaning that somewhere within the structure the contraction of the structure 
slows down due to the pressure like a classical fluid. Therefore, the pressure effect begins somewhere within the 
structure to act like a repulsive force in contrast to the outer regions where its attractive nature dominates. Note 
that the central black hole and its horizon has a much smaller radius than the region of the repulsive pressure effect 
we are discussing. 
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(a)The overall scheme: R > 0 and R < 0 show the expanding and (b)The overall scheme of the acceleration, 

collapsing regions, respectively. 


FIG. 4: The behavior of the Hubble parameter R/R and the acceleration R/R in the case of p = wp. 


C. The Apparent and Quasi Event Horizon 


The boundary of a dynamical black hole, where the area law and the black hole temperature are defined, is a 
non-trivial concept (see for example i, 0 and 13). Our model is again a concrete example to see the behavior of 
both apparent and quasi event horizon of a dynamical structure within an expanding universe. It is easily seen that 
the apparent horizon of our cosmological black hole is located at R = 2M [ij. We note that after the formation 
of an apparent horizon the central density goes to infinity indicating the formation of the singularity at the center. 
Therefore, the singularity is covered by the dynamical horizon and is not naked! 

This apparent horizon is calculated in t, r coordinates numerically. It is always space-like tending to be light-like 
at late times. This can best be seen by comparing the slope of the apparent horizon relative to the light cone at 
every coordinate point of it. This is in contrast to the Schwarzschild black hole horizon where it is always light-like. 
At the late times, however, we expect the apparent horizon to become approximately light-like and approaching the 
event horizon. This is reflected in the Fig.(j5j). It is evident that ^\ah < ^\nuii at all times on the apparent horizon, 
the difference tending to zero at late times. Therefore, the apparent horizon is always a space-like dynamical horizon 
leading to a slowly varying horizon at late times [l,[Il|. Note that the qualitative result is independent of the equation 
of state. 

We now show how the dynamical horizon of our cosmological black hole becomes a slowly evolving horizon at late 
times. Let’s first define the evolution parameter c such that the tangent vector to the dynamical horizon, V, is given 
by 


yn (19) 

where the two vectors and are normal null vectors on a space-like two surface S in (t, r) plane (see 0)- We 
expect c to go to zero at late times in order for our dynamical horizon to become a slowly evolving horizon. In the 
case of our Lemaitre model c is calculated to be 


^ M' + wM' . 
^ ~ ^M' -wM' - 


( 20 ) 


The result of the numerical calculation for different equation of states and different state functions is given in Fig.®. 
The decreasing behavior of the function c in the course of time independent of the equation of state is evident. We 
may then conclude that the dynamical horizon of the cosmological black hole tends to a slowly evolving horizon. 


D. The Effect of A on the Formation of the Apparent Horizon 


Let us now explore how the addition of a cosmological term A to the Einstein equations may influence the formation 
of the apparent horizon. The location of the apparent horizon is now defined by [I^ 

2M = R- . 


( 21 ) 
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(a)The p = wp case: ^\ah < ^\null on the apparent horizon. 
Therefore, the apparent horizon is always a space-like dynamical 
horizon leading to a slowly varying horizon at late times. 



(b)The p = ws{r)p case: ^\ah < ^\null on the apparent horizon. 
Qualitatively, there is no difference to the Fig.(a). 


FIG. 5: 


^\ah < ^\nuii on the apparent horizon. Qualitatively, there is no difference between two equation of states. 



(a)The p = wp case: the more the pressure the sooner the dynamical 
horizon becomes a slowly evolving horizon. 



(b)The p = wps(r) case: qualitatively, the same behavior as in 
Fig.(a). 


FIG. 6: Dynamical horizon becomes a slowly evolving horizon. 


Now, for a typical structure like a galaxy or a cluster of galaxies we have Therefore, the location of the 

apparent horizon is almost the same as in the case of vanishing cosmological constant, i.e. 2M « R. We may then 
refer to Fig ([7)) for the behavior of the apparent horizon. 


E. Mass and matter flux 

Due to the expanding background we expect the matter flux into the dynamical black hole to be decreasing and 
the dynamical horizon to become a slowly evolving horizon in the course of timejl^. We know already that there is 
no unique concept of mass in general relativity corresponding to the Newtonian concept. The question of what does 
general relativity tell us about the mass of a cosmological structure in a dynamical setting was discussed recently 0- 
It was shown [l^ that The Misner-Sharp quasi-local mass, M, is very close to the Newtonian mass. 

Let us then take the Misner-Sharp mass for this black hole and calculate the corresponding matter flux into the 
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(a)The p = wp case. (b)The p = ws{r)p case. 


FIG. 7: Adding A to the Einstein equations does not have any observable effect on the formation of the apparent horizon. 




(a)The p = wp case: the rate of matter flux into the black hole (b)The p = ws{r)p case: qualitatively the same behavior as in 

decreases with the pressure. Fig.(a). 

FIG. 8: The rate of matter flux into the cosmological black hole. 


black hole. In the case of Lemaitremodel, the matter flux is given by 


dM{r,t) dM(r,t) dM{r,t)dr ■ ,dr 


( 22 ) 


The result of the numerical calculation is depicted in Fig.(jH]). Note how the pressure decreases the rate of matter 
flux into the black hole at the late time. 


IV. MEASURING THE REDSHIFT IN THE LEMAITRE METRIC 

In cosmology we are used to interpret the cosmological redshift according to the homogeneous FRW model. What 
if the universe is inhomogeneous? In the simplest case we are ready to model a source, a cosmological structure, 
within an otherwise homogeneous FRW model using our Lemaitre model. We assume now an observer far from the 
source of light near the structure. What is then the redshift measured by this observer? Given that the metric is an 
exact solution of the Einstein equations, we expect the redshift to include all gravitational effects including not only 
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the cosmological redshift but also the gravitational redshift due to the overdensity of the source. 

The redshift in our model can be obtained as follows. Assume the light ray coming from a source S located near to 
the structure in the center of our inhomogeneous Lemaitre model, and the observer O somewhere within the FRW 
background having the corresponding 4-velocities and Let = dx^/dP be the tangent vector to the null 

geodesic connecting the source to the observer. The corresponding redshift z, i.e. frequency shift, is then defined 


as 


22 | 




where and are evaluated at the source and observer events. Let us assume the null geodesics to 

be radial, i.e. = ^ = . For the metric in equation (1) we then have 



62 


dr 


Using now the geodesic equation, it is straightforward to show that 


dk^ 

dp 



where P is an affine parameter. Therefore, 


(24) 


(25) 


dk* 


-[(T' + (-(T+l)^]ci/3fc’’ 


Integrating the Ea. (pHll we obtain 





Co exp (- 



(26) 


(27) 


where Cq is a constant. 

Now, the 4-velocities are given by 



(28) 

uy = 

(29) 


Using equations (IMl) . (l28ll . (IMl) . and (1271) we obtain the cosmological redshift in the presence of a structure {zcbh)' 

1 + zcBH = exp[{j^ + i-<^+^)^]dr)]^. (30) 


1 + ZCBH = eMif [(-^ + ^)p]dr)] = eMif ip/p)]dr)]. (31) 

Note that there was no need to calculate the k'^ due to k^u^ = 0. 

Now, the Eq. (| 3T|) may be integrated numericaly using the equations for d, cr', and A from [2^ for any collapsing 
structure in an expanding FRW universe. The necessary initial condition may be chosen along the line discussed in 

0. 

This redshift includes the familiar cosmological FRW part, zc, as well as the gravitational redshift, zq due to the 
overdensity of the structure (the cosmological black hole). In general, we then expect it to be different from that of 
the corresponding homogeneous FRW model. In the special case of a homogeneous universe without a structure it 
reduces obviously to the familiar FRW cosmological redshift zc = — 1, lacking the contribution from the local 

CiyZs ) 
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gravitational redshift of the overdensity, zq- To see the difference between the exact inhomogeneous cosmological 
redshift according to our model and the sum zq let’s look at some specific models. 

Take a CBH model with the mass 10 ^Mq at a distance corresponding to zc = 0.005 from the observer. The 
gravitational redshift according to the Schwarzschild metric is given by 




1 


( 1 - 


2M(r,t) 

Rs 


1 . 


Adding to it the FRW cosmological redshift 




a{to) 

aits) 


should give us the CBH redshift we have already calculated, i.e. 


(32) 


(33) 


ZCBH Z:^Zg + Zc (34) 

The result of the numerical calculation is given in Fig. We see that the exact zcBH is always larger that the 
sum of FRW cosmological redshift and the local Schwarzschild gravitational redshift, although the difference is smaller 
than the observational limit of accuracy. The difference zcBH — zc is shown in Fig. (pR|) . The difference goes to zero 
for the source being at distances far from the apparent horizon. There may be, however, cases that this difference is 
not to be ignored. We leave it to a more detail study in future to see where this difference may be of any cosmological 
significance. 



FIG. 9: Difference of gravitational redshift and CBH redshift for the some source near the CBH {Rah — Rapparent-horizon)- 
The observer place is fixed. 


R/Rah 

ZG ZCBH [zCBH 

— {zG + Zc)]/{ zcBh) 

30 

0.01709 0.022173 

0.00374 

50 

0.01015 0.015211 

0.00401 

70 

0.00722 0.012263 

0.00350 

100 

0.00458 0.00959 

0.00101 


TABLE I: Numerical values of Fig.([^ 


V. MODELING THE EVOLUTION OF A CLUSTER OF GALAXIES USING THE LEMAITRE 

METRIC: CASE OF A2061 


Now we are ready to use our exact solution and the corresponding algorithm to study the collapse of a real galaxy 
cluster within an otherwise expanding universe using a fluid model. Clusters of galaxies are modeled by a spherically 
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FIG. 10: The difference between cosmological and CBH redshift for several observers located at different physical distances, 
the source of photon is fixed for all observers. 


symmetric (isotropic) dark matter halo assumed to dominate the dynamics of the system. The density profile of such 
a halo is often described by the Navarro-Frenk-White (NFW) profile M- 


PNFw{r) 


^cPc 


^(1+^)2’ 


(35) 


where Vg is a scale radius, pc is the critical density, and the characteristic over density 5 c being a function of cluster 
concentration c = r2oo /given by 


200 

3 ln(l + c) - ■ 


(36) 


The parameter r 2 oo is the radius at which the average interior density is 200 pc being approximately equal to the 
virialized over density, pviriai = 178 pc- The NFW profile in Eq.(j35j) has three undesirable features to be applied to 
our case of a cosmological structure: 

1) The density tends to infinity at small r. 

2) The NFW mass diverges at large r and can not be matched to the with FRW mass; 

3) Due to the lack of a void it doesn’t represent a density profile for a structure within an expanding universe. 

We have then to modify it to remedy these deficiencies. To resolve the first problem, we introduce a maximum density 
at very small r: 


pi NFw{r) = 


^cPc 


(ec + ^)(l + ^)2- 

The second problem will be resolved by introducing the truncation radius r^: 


p2 NFw{'^) 


5cPc 


(ec + ^)(1 


K: 


.)2 '•^2 


(37) 


(38) 


In order to resolve the third problem, we add a Gaussian density profile to match the overdensity remon through a 
void to the cosmological background as is necessary for any exact solution of Einstein equations (see [6| and 0 ): 


pair) = aexp {-——^^). (39) 

To 

The parameters ro and ri give us the freedom to adjust the location of the void as desired. Eor our model cluster, 
A2061 [26|, we choose c = 3, a = IMpc, and ec = 0.1 to have r 2 oo = 1.723Mpc . Our metric solution has to approach 
ERW at large r outside the cluster. Therefore, we add a homogeneous background density such that the final density 
profile at to becomes 


Pcluster{r) = p2 NFw{r) — Pg{t) + pc- 

Note that the Gaussian profile must be added in such a way that 


pri on 

/ r‘^pQdr= / r^p2{r)dr, 

Jo Jo 


(40) 


( 41 ) 
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where ri is the distance at which the density approaches the background critical density. It guaranties that the mass 
due to the overdensity in the structure is compensated by the underdensity region within the void (See Fig. (pT]) b 
Now, by choosing i?(to,r) = r and integrating Eq. ^ to determine M(to,r) as well as having the equation of state, 
we can solve the coupled evolution PDEs numerically. To have zero pressure at ERW background, a good choice for 
the equation of state is p(t, r) = w{p{t, r) — pc), in which pc is just a function of time. The other initial functions are 
as before. Eig. m shows the evolution of the galaxy cluster. By modifying the NEW density profile to match the 
exact general relativistic requirements we have therefore arrived at a model cluster based on the exact solution. We 
may even go further and see how good we may rely on the observation of the red shift of such a cluster and compare 
it to the model reflected in Eq. (|3T]) . 

Based on the redshift of clusters estimated with spectroscopic data and the redshift of the brightest cluster galaxies 
(BCGs) [26| we may take the cosmological redshift of the galaxy cluster A2061 to be 2: = 0.07845, being located 
at a distance around 330 Mpc. We may therefore model our galaxy cluster such that its cosmological redshift is 
= 0.07845; The Eq. (|3T]) will then tell us if we need to add the gravitational redshift as part of the cosmological 
one. 

Eixing the observer distance at 330Mpc, and solving the Eq. (|3T]) we see that the gravitational redshift is about two 
orders of magnitude smaller than the cosmological redshift. Therefore, we have a well defined dynamical model of 
the structure including a modified density profile fitted well with observation. The effect of gravitational redshift 
may be much bigger for supermassive blackholes [25| . Any other density profile may also be modified and adapted to 
requirements as indicated in this section. 



FIG. 11: Density evolution of the galaxy cluster. 


VI. DISCUSSION 

We have studied the evolution of a structure made of perfect fluid with non-vanishing pressure as an exact solution 
of Einstein equations within an otherwise expanding ERW universe. The structure boundary is separated by a void 
from the expanding part of the model which is very much like a ERW universe already near by the void. We have 
noticed a counter-intuitive pressure effect somewhere inside the structure where the existence of the pressure slows 
down the collapse like a classical fluid in contrast to distances far from the structure. The collapsed region develops to 
a dynamical black hole with a space-like apparent horizon, in contrast to the Schwarzschild black hole. This apparent 
horizon tends to a slowly evolving horizon becoming light-like at late times with a decreasing mater flux into the 
black hole. We have, therefore, to conclude that the mere existence of a cosmological matter, even dust, may have 
significant effect on the central black hole differentiating it from a Schwarzschild one irrespective of how small the 
density outside the structure is. 

The light properties of these cosmological black holes can be interesting because most of the information about the 
black holes, galaxies and clusters are coming from their lights. We have investigated the redshift of a light emitted 
near by a cosmological structure to a distant observer. It was shown that the exact CBH redshift of light includes 
both local gravitational and cosmological redshift of the structure. Therefore, in the era of precision cosmology we 
may be forced to consider the effects of inhomogeneities as seen in the CBH cosmological model for high redshift 
surveys. 

We have also seen how to generalize the existing models of density profile for cold dark matter within large structures 
using the results of our structure model. Although the gravity may be too weak near a large cosmological structure. 
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we can not use the Newtonian approximation due to the non-local or quasi-local cosmological effects. 
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